Effects of Interactions and Temperature in Disordered Ultra-Cold Bose Gases 
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We simulate ultra-cold interacting Bosons in quasi-one-dimensional, incommensurate optical lat- 
tices. In the tight-binding limit, these lattices have pseudo-random on-site energies and thus can 
potentially lead to Anderson localization. We explore the parameter regimes that lead to Ander- 
son localization and investigate the role of repulsive interactions, harmonic confinement and finite 
temperature. We find that interactions can obscure the exponential localization characteristic of 
Anderson localization, thus impeding the direct observation of this phenomenon when interactions 
are present. 
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Atoms in optical lattices provide a great tool to study 
many theoretical models from solid-state physics as they 
allow for unique control of normally unmodifiable param- 
eters. Disorder can be introduced and controlled by ex- 
ternal optical potentials [l|, and interactions can be 
tuned by a Feshbach resonance Q ■ While being intrinsic 
to a solid-state system, disorder can easily be imposed 
onto and modified in an optical lattice. Disorder in a 
quantum mechanical system leads to various interesting 
effects such as Anderson localization, which describes the 
localization of matter waves in a random potential due to 
destructive interference from randomly scattered waves 
[J]. Several experimental groups have studied localiza- 
tion effects in ultra-cold Bose gases in disordered optical 
lattices. To create the disorder, either a laser speckle 
potential Q or a second optical lattice incommensurate 
with the original lattice is superimposed onto it d,0|. We 
extend mean field || treatments and their finite temper- 
ature generalizations 0-[H| > which have been extensively 
developed for the dilute Bose gas, including their appli- 
cation to the optical lattice [12|, to include disorder. 

We study a Bose gas in a quasi-one-dimcnsional disor- 
dered optical lattice. A tight-binding model leads to a 
Bose-Hubbard Hamiltonian fl2| which is further approx- 
imated and numerically solved. We study the effects of 
different disorder and interaction strengths on the con- 
densate density, the condensate and superfluid fractions, 
the lowest excitation energy and the band structure. We 
use the same setup as Schulte et al. [7$, but also treat a 
harmonically confined Bose gas at finite temperature. 



I. MODEL 

A. The system 

The model we use for our simulations is taken from 
[l2| and references therein. Tight radial confinement of a 
Bose gas leads to a quasi-one dimensional problem while 



scattering remains three-dimensional. The lattice poten- 
tial is Vl{x) — Vqcos 2 {2ttx/\l)i where Vq is the lattice 
depth and Al is the lattice wavelength. The recoil en- 
ergy is defined as En = f^ki 2 /2m, where kh is the lattice 
wave-number. The disorder potential is of the form 
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where a and j3 are the incommensurate ratios between 
the primary and secondary lattice wavelengths. For irra- 
tional values a and f3, the resulting potential is pseudo- 
random. In experiment, however, the wavelength ratios 
arc always rational, so the potential exhibits an addi- 
tional level of quasi-periodicity. If a and /3 are ratios of 
large coprime integers, however, the period of the poten- 
tial will be larger than the system size. 



B. The Bose-Hubbard Model 

The wavefunction is expanded in Wannicr functions lo- 
calized at each lattice point of the lattice potential. This 
tight-binding approximation leads to the Bose-Hubbard 
Hamiltonian for atoms in a one-dimensional optical lat- 
tice of N sites, 
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with tunneling energies J; i? :+i, on-site interaction poten- 
tial U, and on-site energy e^, which contains both the 
trapping potential and the disorder potential. These 
parameters are determined by the lattice and disorder 
potential and are obtained by taking integrals over the 
Wannier functions [l6j . The resulting expressions in the 
Gaussian approximation [T3] are 
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a s is the scattering length of the 



atoms and coj_ the radial (perpendicular) trapping fre- 
quency. 

The disorder potential leads to pseudo-random shifts in 
the on-site energies and the tunneling energies. The on- 
site energy for one secondary lattice with wavelength ra- 
tio a to the primary lattice wavelength is 



A. Homogeneous Trap at Zero Temperature 

To isolate the effects of disorder and interactions, the 
first simulations are performed without a harmonic trap- 
ping potential. We refer to this as the homogeneous case 
albeit dealing with a finite size system with hard wall 
boundary conditions. Fig. [1] and [2] show the condensate 
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for two secondary lattices with ratios a and /3, the on- 
site energy is just the sum of the above term for each 
of the two ratios. Our simulations show that the site- 
dependence of the tunneling rate has a negligible effect 
on the condensate and we thus omit it here for simplicity, 
such that Jis+i = J for any i. 



C. The Hartree-Fock-Bogoliubov Formalism 

When a macroscopic occupation of the ground state is 
assumed, the Bose annihilation operator can be written 
as a, = (zi + 5i)e~' l ^ t / h , where z, is a complex mean-field 
part and Si the fluctuation operator [l2j . The condensate 
density is given by n c ,i = \zi\ 2 and the non-condensate 
density by n t = (SjSi). 

From this, we can construct a discretized, generalized 
Gross-Pitaevskii equation 



[izi = tiz t - J (zi+i + Zi-i) + U(n c ,i + 2hiZi), 
and Bogoliubov-de Gennes equations 
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and v\ are the Bogoliubov quasi-particle amplitudes and 
Nbe is the Bose- Einstein distribution. This is a lat- 



tice formulation [12| 
Fock-Bogoliubov [! 



13i of the finite-temperature Hartrec- 
4llj formalism. The equations are 



then solved sclf-consistently. 



II. RESULTS 

Here we show results for total atom number N = 500 
atoms of 87 Rb in a system with 100 lattice sites in the 
experimental configuration of 0, where a = 55/64 and 
(3 = 165/212. 

The strength of inter-particle interactions can be tuned 
by a Feshbach resonance Q and, for finer control, by the 
radial confinement. Here we define A = eo/J and V e ff = 
U j J so that we need not be concerned with evaluating 
the integrals from experimental parameters. 
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FIG. 1: Condensate density for various parameter 
values, (a-c) V eff = 10~ 8 . (d) V eff =6x 10" 4 . 



10° 



io-' 



10° 



s 

■o 10 -io. 



-40 

(a) A = 





site 



40 



-40 




site 



40 



(b) A = 0.68. 




(c) A = 1.13. 



(d) A = 1.13. 



FIG. 2: Condensate density (logarithmic) for various 
parameter values, 
(a-c) V efJ = IO" 8 . (d) V e ff = 6 x 10- 4 . 

density for different parameters in linear and logarith- 
mic scale. The first panel shows the weakly interacting, 
non-disordered case, which is just the ground-state of a 
single particle in a box. Increased disorder then leads 
to a narrowly peaked density. The logarithmic plot in 
Fig. He] indicates an exponential decay of the conden- 
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sate density, a signature of Anderson localization. The 
Andre- Aubry model [bij , which describes non-interacting 
Bosons in two incommensurate lattices, predicts a transi- 
tion from extended to localized states for A = 2. Here we 
have two secondary lattices instead of one, so the tran- 
sition occurs for A = 1. In the presence of interactions, 
the condensate fragments into multiple peaks. Each sin- 
gle peak remains exponentially localized (Fig. I2d|) . but 
limited spatial resolution in experiments might be unable 
to distinguish them. 

We find that increasing the disorder strength contin- 
uously lowers the superfluid fraction, underlining the lo- 
calizing effect of disorder (Fig. [3]). The decrease is not 
dramatic as even for strong disorder the superfluid frac- 
tion remains above 0.85. The long-range phase coherence 
implied by this may well be an artifact of our mean field 
technique, which is not strictly valid in the localized ex- 
treme. The technique should, however, be reliable up to 
this point and is thus a valid predictor of any transition. 
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FIG. 3: Superfluid fraction f s vs. A at various 
interaction strengths. 

In the case of weak disorder, the superfluid fraction 
decreases slightly with increased interaction strength, 
which is explained by repulsive interactions inhibiting 
superfluidity. In contrast, for increased disorder the su- 
peruid fraction increases with higher V e ff ■ This fol- 
lows from the interplay between interactions and disor- 
der. The disorder localizes the condensate, thereby re- 
ducing the superuid fraction. The interactions delocalize 
the condensate and thus attenuate this effect. We note 
that for very large interaction strengths, the superfluid 
fraction will decrease again [TBI ]. 

Both the superfluid fraction and the lowest excitation 
energy show a sharp drop in the region where the in- 
teractions begin to fragment the condensate (Fig. @|. 
The lowest excitation energy approaches zero where the 
energies for the single-peaked state and the fragmented 
state are equal. This effect results from the competition 
between the on-site energy and the kinetic and interac- 
tion energies. The single large peak minimizes the ex- 
ternal potential energy while a fragmented condensate 
minimizes the kinetic and interaction energy Above a 
critical interaction strength, the increased potential en- 




FIG. 4: Superfluid fraction f s and lowest excitation 
energy E\ vs. U at A = 1.13. 



crgy is compensated by the reduced kinetic and interac- 
tion energies. As each peak is exponentially localized, 
the superfluid fraction decreases at the critical interac- 
tion strength, but afterwards continues to increase as the 
interactions keep delocalizing the condensate. 

The excitation spectrum is shown in Fig. [S]for a weakly 
interacting condensate with and without disorder. In the 
absence of disorder, we obtain the band structure of a 
free particle in a sinusoidal potential. In the presence of 
disorder, the band shows gaps at fractions corresponding 
to the lattice ratios a and (3, demonstrating the quasi- 
periodicity of disorder. In the setup of 0, for example, 
a is approximately 0.86 and /? is approximately 0.78. We 
find gaps, for example, at indices 86 and 100 — 86 = 14 
due to a, at 78 and 100 — 78 = 22 due to (3 and at 35 
due to a beating, where 0.35 = 2(a~ 1 + (1 — 
Our simulations show that V e f / has no qualitative effect 
on the band structure. We conclude that the excitation 
spectrum is primarily governed by single-particle effects. 




FIG. 5: Band structure of the weakly interacting 
condensate (V e ff = 10~ 8 ). 
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B. Harmonic Trap 



C. Finite Temperature 



The general observations of the previous section re- 
main valid when a harmonic trapping potential is present, 
with the main difference being the additional confinement 
due to the trap. In this case, only in the logarithmic plots 
is it possible to distinguish between localization due to 
the trap, which is Gaussian, and localization due to An- 
derson localization, which is exponential (Fig. |6]). Wc 
also note that the trap inhibits the fragmentation of the 
condensate. This can be seen most clearly in the graphs 
for the superfiuid fraction and the lowest excitation en- 
ergy (Fig. [7]). In contrast to the homogeneous case, no 
sharp drops occur and the lowest excitation energy does 
not approach zero. Wc conclude that no sudden crossing 
from a single-peak state to a multiple-peak state occurs. 
Instead, side peaks arise continuously over a wide range 
of interaction strengths. 




(c) A = 1.13 



(d) A = 1.13 



FIG. 6: Condensate density (logarithmic) at f t = 10Hz 
for various parameters, (solid black) V e ff = 10 -8 . 
(dashed red) V eff = 3.7 x 1CT 4 . 
(dotted blue) V e ff = 6.2 x 1CT 4 . 
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FIG. 7: Superfiuid fraction and lowest excitation energy 
vs. interaction strength. 



We now discuss the situation for finite temperature. 
For the mean-field approximation to remain valid, wc 
must stay well below the critical temperature where the 
non-condensate density is low. For our simulations, we 
take the setup from the previous section and vary the 
temperature to lOnK. 

The effect of the finite temperature is to drive atoms 
out of the condensate into the thermal cloud. Thus the 
condensate fraction f c is lowered and the non-condensate 
density h is increased. Because Anderson localization is 
a single-particle effect that docs not depend on the to- 
tal number of atoms in the condensate, we still see a 
Gaussian peak for weak disorder and exponential local- 
ization for strong disorder. The thermal cloud, while be- 
ing slightly affected by the disorder potential, does not 
become localized. This is because the thermal cloud is 
incoherent and therefore behaves more like a classical gas. 
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FIG. 8: Condensate fraction vs. interaction strength at 
T = lOnK. 

The interplay between disorder and interactions is 
more subtle at finite temperature. We first observe that 
increasing the interaction strength does not fragment the 
condensate. Instead, atoms are driven into the thermal 
cloud (Fig. [5]). This is energetically advantageous as the 
thermal cloud has lower interaction energy due to its re- 
duced density. We further see that disorder increases the 
condensate fraction for weak interactions (Fig. [5]) and de- 
creases it for strong interactions. The former is a conse- 
quence of micro-trapping. The disordered potential con- 
tains microscopic wells whose harmonic trapping frequen- 
cies are higher than the 10Hz of the external trap. This 
increases the critical temperature Tq, because Tq ~ 
for a BEC in a harmonic trap of trapping frequency uj t - 
Consequently, the ratio T/Tq is decreased and the con- 
densate fraction rises. For strong interactions, in con- 
trast, the localizing effect of disorder increases the in- 
teraction energy, thereby driving more atoms into the 
thermal cloud. 

The superfiuid fraction follows the condensate fraction 
because only atoms in the condensate contribute to the 
superfluidity. Apart from this, it does not show any new 
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behavior differing from that of the zero-temperature case. 
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FIG. 9: Condensate fraction vs. disorder strength at 
T = lOnK and V eff = 10~ 8 . 



D. Excitation Energies 

We conclude the treatment of the harmonically con- 
fined Bose gas at finite temperature with a discussion of 
the excitation energies. Again, we look at the interplay 
of temperature, disorder and interactions. We expect 
the behavior of E\ to correspond to the behavior of f c . 
Lower excitation energies should result in a lower con- 
densate fraction because less thermal energy is needed to 
populate low lying levels. 

First, we investigate the temperature dependence. Our 
simulations show that for weak interactions (V e ff ~ 
10 -8 ), the temperature has essentially no effect on the 
band structure, regardless of the value of A. This is not 
surprising as temperature does not alter the energy lev- 
els or states of a single particle system. In the presence 
of interactions, however, we see that increasing the tem- 
perature also increases the lowest excitation energy (Fig. 
[TO]) . and this effect becomes more pronounced as V e ff 
increases. This can be explained with the simultaneous 
decrease in the condensate fraction as temperature and 
interaction strength increase. Removing atoms from the 
condensate lowers the energy of the ground state via the 
interaction contribution n c V e ff and raises the energy of 
the excited states, thereby increasing the energy differ- 
ence Ei between the ground state and the first excited 
state. It is worth noting that, contrary to our expec- 
tation, the condensate fraction does not increase with 
increasing excitation energy. This is attributed to the 
dependence of E\ on f c in the presence of interactions, 
as discussed above. 

Next, we turn to the interplay between disorder and 
interactions. Fig. [TT]depicts the lowest excitation energy 
vs. V e ff. The graph shows that in the presence of dis- 
order, increased interactions lower the excitation energy. 
At constant temperature, this allows more atoms to 
leave the condensate, in agreement with what we observe 
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FIG. 10: Ei vs. T for various interaction strengths at 
A = 0. 
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FIG. 11: Ei vs. U for various disorder strengths at 
T = 10 nK. 



for f c . The effect of disorder in the weakly interacting 
regime is to increase the excitation energy. For higher 
interaction strengths, this effect is reversed. This reflects 
the fact that more atoms leave the condensate when it 
becomes localized in the presence of strong interactions. 



III. CONCLUSION 

In this paper, we have discussed the interacting ultra- 
cold Bose gas in incommensurate optical lattices and 
highlighted the important role of disorder, interactions 
and temperature in the quest for Anderson localization. 
Our simulations verify that Anderson localization occurs 
and elucidate how finite temperature, a harmonic trap- 
ping potential and repulsive interactions affect the con- 
densate. 



IV. ACKNOWLEDGMENTS 

This work was supported by contract SFB/TR 12 
of the German Research Foundation and through the 



6 



IB BMBF (Project NZL 07/006) and by the New nologies and the New Zealand International Science and 
Zealand Foundation for Research, Science and Technol- Technology Linkages Fund, 
ogy through contract NERF-UOOX0703:Quantum Tech- 



[1] M. Lewenstein, et al, Adv. Phys. 56, 243 (2007). 
[2] L. Fallani, et al, Adv. Atom. Mol. Opt. Phy. 56, 119 
(2008). 

[3] J. M. Vogels et al, Phys. Rev. A 56, R1067 (1997). 

[4] P. W. Anderson, Phys. Rev. 109, 1492 (1958); for a re- 
view, see; B. Kramer and A. MacKinnon, Rep. Prog. 
Phys. 56, 1469 (1993). 

[5] L. Sanchez-Palencia, et al, Phys. Rev. Lett. 98, 210401 
(2007). 

[6] L. Fallani, et al, Phys. Rev. Lett. 98, 130404 (2007). 

[7] T. Schulte, et al, New J. Phys. 8, 230 (2006). 

[8] C. J. Pethick and H. Smith, Bose-Einstein Condensation 

in Dilute Gases, Cambridge University Press (2002). 
[9] A. Griffin, Phys. Rev. B 53, 9341 (1996). 
[10] D. A. W. Hutchinson, E. Zaremba and A. Griffin, Phys. 

Rev. Lett. 78, 1842 (1997). 
[11] D. A. W. Hutchinson, et al., J. Phys. B 33, 3825 (2000). 



[12] B. G. Wild, P. B. Blakie and D. A. W. Hutchinson, Phys. 

Rev. A 73, 023604 (2006). 
[13] A. M. Rey, et al, J. Phys. B 36, 825 (2003). 
[14] S. Aubry and G. Andre, Ann. Israel Phys. Soc. 3, 133 

(1980). 

[15] G. G. Batrouni and R. T. Scallettar Phys. Rev. B 46, 
9051 (1992). 

[16] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 
(2006). 

[17] In this approximation, the Wannier functions in the 
expectation-value integrals are replaced by Gaussians. 
This gives accurate results for the on-site energies while 
underestimating the tunneling energy. Since results are 
given in terms of the tunneling energy J, this is of no 
concern if one keeps in mind that a given value for J does 
not exactly correspond to the given lattice geometry. 



